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Abstract 

The purpose of this article is to study the persistence of solution of a 
hyperbolic system under small viscous perturbation. Here, the solution of the 
hyperbolic system is supposed to be periodic: it is a periodic perturbation 
of a roll-wave. So, it has an infinity of shocks. The proof of the persistence 
is based on an expansion of the viscous solution and estimates on Green's 
functions. 
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1 Introduction 

In this paper, we consider a one-dimensional system 



e 

XX 



with a smooth flux / : R" — )■ R" and a smooth function g : R" — )■ M". We assume 
that the corresponding system without viscosity 

ut + fiu)x = giu) (2) 

is strictly hyperbolic. 

We consider a piecewise smooth function u which is a distributional solution 
of ([2]) on the domain R x [0;T*]. We assume that u is periodic in the x variable, 
with a period L and that u has m noninteracting Lax shocks per period. 
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Introduction 



We show that m is a strong hmit of solutions of ([T]) as £ — ?■ 0. This work is of 
course motivated by the conjecture that the admissible solutions of ([2]) are strong 
limits of solutions of ([T]) with the same initial data. 

In the case of scalar conservation laws, the proof of this conjecture uses the 
maximum principle ^U\, and in the case of special 2x2 systems, R. J. DiPerna proved 
it by a compensated compactness argument [1]. For the general case of shocks, 
there is a first paper of J. Goodman and Z. P. Xin which proves this conjecture 
for small amplitude Lax shocks [3]. This conjecture is also proved for a single non- 
characteristic Lax shock or overcompressive shock by F. Rousset [9j. Here, we only 
consider Lax shocks but we have an infinity of shocks. 

An other motivation of this work states in the study of roll-waves, in fluid me- 
chanics or in general hyperbolic systems with source terms. Indeed, P. Noble proved 
the existence of roll-waves for this kind of system under assumptions on the source 
term [?]. Specifically, in the case of inviscid Saint Venant equations 



one can prove that there exist roll-waves which are persistent under small perturba- 
tion [8]. So, there exist solutions of inviscid Saint Venant equations, near roll- waves. 
Here, the idea is to prove that there exists a family of solutions of the viscous Saint 
Venant system 



which tends to a solution of inviscid system as e goes to 0. We prove this result in 
the case of full viscosity. 

We can now give the full set of assumptions and formulate our main result. First, 
we suppose that 

(HI) system ([2]) is strictly hyperbolic. 

That means that there exist smooth matrices P{u), D{u) such that 



where D{u) = diag(Ai(M), . . . , A„(-u)) is a diagonal matrix and Aj 7^ \j for all i 7^ j. 

(H2) M is a distributional solution of on [0; T*]. Moreover, we suppose that u is 
piecewise smooth, periodic, and has m noninteracting and non-characteristic 
Lax shocks per period. 





(3) 



df{u) = P{u)D{u)P{u)-^ 



2 




m m 711 7n 

Xi{t) X2{t) X^{t) Xi{t) Xi{t)+L 



Figure 1: Allure of solution u over one period when u is scalar and m = 4. The 
periodic roll-wave is drawn in dotted line. The solution which checks our assumptions 
is represented by continuous line. One also placed the shocks for the two solutions. 

That means that u is smooth except at the points {x, t) of smooth curves x = 
Xj{t) + iL,j = 1, . . . ,m,i E Z and that for all j, k, t, \Xj{t) — Xk{t)\ > 2r > (see 
Figure [I]) . 

Moreover, following limits are finite: 

ay ±(t) := d'MX,{t) ±0,t)= hm d'M^, t). 
Since the shocks are non-characteristic k-Lax shocks, we have: 

Xi{u'-) <■■■< Xk-i{un < x'^it) < <■■■< K{un, 

Ai(m^+) < ■ ■ ■ < \k{u^+) < x;(t) < Afc+i(M^+) < ■ • • < A„(m^+). 

This assumption ensures the existence of at least one sonic point between two shocks. 

We refer to [8] for the existence of such a solution in the case of Saint Venant 
equations. This result can be extended to general hyperbolic systems. 

These assumptions imply that there exists a viscous profile for each shock. More 
precisely, for all j, there exists such that 

v^^-{fiv^)-x;v^)^ = o (4) 

and 

V^{±oo,t) = u^^{t). 

We will give more details on the properties of in Section !!^. 2.11 Now, we only need 
to expose some assumption of linear stability. Consider for t <T*, the operator 

£> = w,, - {df{V^{z, r)) - X'^{t))w,. 
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We assume that the viscous shock profiles are hnearly stable. This assumption is 
equivalent to an Evans function criterion [TT] . 

(H3) Vr e [0; T*], j = 1, . . . , m, Ci is such that D^(A) ^0 VA, 3ftA > 0, A 7^ 0, and 
Dl (0) 7^ 0, where is the Evans function of Cl. 

We can now state our main theorem: 

Theorem 1. Under assumptions (HI)— (H2)— (H3), for all e > 0, there exists a 
unique solution of ([T]) on [0; T*] such that 

u'{t = 0,x) =u{t = 0,x). (5) 

Moreover, we have the convergences 

- ^^||l-'([0;T*],L1(0;L)) -^0, aS E 0. 

And for any rj G (0, 1), 

sup \u'^{x,t) — u{x,t)\ ^ 0, as £ — )■ 0. 

0<t<T'-,\x-Xj(t)\>e^ 

The proof of this theorem is done in three steps: construction of an approximate 
solution (which gives us an expansion of u'^ in e), estimates on the semigroup gen- 
erated by linearized operator around this approximate solution and a Banach fixed 
point argument to deal with the full problem. 

The paper is organized as follows. In Section O we build an approximate solution 
u^pp of the full problem ([T]) close to u, solution of ([2]) up to order 2 with respect to 
e. This is done separating slow parts where m^^^ is close to u and shock parts where 
ulpp\xj±e-' is close to Vj. More precisely, one expands in the slow part as 

t) = u{x, t) + eui{x, t) + e'^U2{x, t) + oi^e^) 

where u is the solution of and Ui are solutions of the linearized equation of ^ 
around u, which is well-posed thanks to assumption (H3). In shock parts, the 
expansion at shock j is 

u'{x,t) = V^{e{x,t,e),t) + eVl{e{x,t,e),t) + e^Vi{e{x,t,e),t) + o{e^) 

where the stretched variable is ^^{x,t,e) = ^ + S^{t), is solution of viscous 
equation (jl]) and V/,i = 1,2, are solutions of the linearized equation of (jl]) around 
. Moreover, the functions Ui, V/ are related by matching conditions which ensure 
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regularity on the approximate solution , built by convex combination of the 



expansions: 



where 



if |a;| > 2, 

1 if Ixl < 1, 



P%x,t) = V\e{x,t,e),t)+eVUe{x,t,e),t) + eX{e{x,t,e),t), 
0^(x, t) = u{x, t) + eui{x, t) + e^U2{x, t). 

With this construction, we prove the theorem: 

Theorem 2. There exists an approximate solution w^^p of ([T]) defined on [0;T*]. // 
if is a smooth change of variable which fixes the shocks (\/t,i,j,{p{{j — l)^ + iL,t) = 
Xj{t) + iL), and ulpp{z,t) = ulpp{(p{z,t),t), then ul^p verifies the equation 

iKpp)t + f{ulpp)x - e{ulpp)xx - giulpp) = t) 
with the following estimates on (f 

ll^t llQ'tilU°° < Ce^^ , (6) 

W\\L^(<d;L)i llg^i ||l1(0;L), II IIli (0;^) ^ Cs'^'^ , (7) 

Il9^lk°°) ll^2tlU°° < Ce'^, Wz\\l^(0]L), 11^2*11x1(0;^) < Ce^'^, (8) 

\WJ\L-<CAWJWio;L)<Ce\ (9) 

Here, w^^p is constructed as a perturbation of u going to order 2, which allows 
us to have estimates on ql^ in in e'^. This property will be useful to prove the 
convergence of m^^^ — to 0. 

In Section [3], we linearize ([I]) in the neighbourhood of the approximate solution 
u%pp and we compute estimates on the Green's function. To do so, we use the method 
of iterative construction of the Green's function, which was first used by E. Grenier 
and F. Rousset in [3]. So, we consider approximations of the Green's functions in 
neighbourhood of the shocks (given by K. Zumbrun and P. Howard in [IT]) and 
we build our own approximation far away from the shock, using the characteristic 
curves. 

The last section is dedicated to the proof of theorem 

Theorem 3. Under assumptions (HI)— (H2)— (H3), for alls, there exists u'^ solu- 
tion of ([I])-© on (0;T*). And this u'^ verifies the convergences: 



W - <pplU°°((0;T*)xM) 0, 
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when e goes to zero. 

This is done using standard arguments for parabolic problems. Indeed, we com- 
bine estimates on g^, and estimate on the Green's function to obtain estimates on 
u^ — M^pp, and its derivatives, depending on e and uniform in time for e small enough. 
Then, using the convergence of u^^p to m, we immediately deduce Theorem [H 

2 Construction of the approximate solution 

The purpose of this section is to prove Theorem|2]on the existence of the approximate 
solution M^pp of ([I]). In a first step, we compute formally this approximate solution 
using outer and inner expansions of order 2. Indeed, in slow part, where Vu is 
bounded, the solution of ([T]) may be approximated by truncation of the formal 
series 

■u^(x, t) ~ 0^(x, t) = u{x, t) + 6Ui{x, t) + e'^U2{x, t) 

where u is the solution of (E]) we want to approach. Similarly, near the shocks j, we 
search for w^pp with the inner expansion 

P^{x,t) = V^{e{x,t,e),t)+eVi{e{x,t,e),t)+e^Vi{e{x,t,e),t) 

where C,-'{x,t,e) = ^'^j'^^^ -)_ 5i{t) + s6{{t) is the stretched variable and is the 
viscous shock profile, solution of (jl]). We match this expansion by continuity of -u^pp 
and its spatial derivatives. 

In this section, we formally substitute these expansions in ^ to find equations 
satisfied by Ui and y/, i = 1,2, j = 1, . . . ,m, and matching conditions. Then, we 
prove the existence of the and V/ on (0; T*). Furthermore, we give rigorous esti- 
mates on the error terms. We can remark here that we search for an approximation 
of order 2, this will be useful to obtain good estimates on the second derivatives of 
the error term. 

2.1 Formal calculation and derivation of the equations 

Substituting into ([T]) and identifying the power of e in the expressions, we get 
for X ^ Xj{ty. 

O(£°):ni + (/(n)).-(7H = 0, 

0{e^) : ui^t + (d/(M) • - dg{u) ■ ui = u^^, 

0{e^) : U2,t + (d/(M) • U2)x - dg{u) ■ U2 = ui^^ - ]^{<l^f{u) ■ {ui, mi))^: 

+ \ '^^9{u) ■ {ui,ui). 
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A similar calculation for yields the set of equations: 
O{e-'):V^^~{f{V^)-X'^V% = 0, 
0{e') : ^4 - miVn - Xj) ■ V^)^ = + V^6i, - g{Vn, 

We remark that the equations for Uj are hyperbolic equations: the first one 
is (E]), so nonlinear, and the others are the linearization of (E]) around u. Similarly, 
the equations for the shock profiles are ordinary equations: nonlinear for we 
recognize (jl]), and its linearization around for Vl and V^. To maximize the 
order of the approximation, we couple these equations with boundary conditions, 
connecting Ui and Vl . First, we note 

dlu\^{t) := dlulX,{t)±^,t) = \im d'Mx,t), t = 1,2. 

Then, we rewrite and P with the variable ^, in a vicinity of shock j, and we ask 
the two functions to coincide as e goes to 0. Therefore, we make Taylor expansion 
of order 2 with respect to e. For example, for > 0, large enough, 

0'iX,{t) + - 5i{t) - s5i{t)),t) = u^^t) +s{u{^it) + ui^i^ - 51)) 

+ y K^i) + 2«i:(t)(^ - - + " ^1?) + o{e') 

and 

r{X,{t) + e{i - m - e6m,t) = V^iU) + eViiU) + e'VUU) + o{e') 
Identifying the terms of same order on e, we get as ^ — > ±oo: 

{±oo, t) = u^^{t), (10) 

t) = utit) + uiHm - m)) + (n) 
VI t) = u^{t) + u^^im - m) + l^itim - mr - <Ht)m + o(i). 

(12) 

For more details on the computation of these conditions, we refer to [2]. 

2.2 Existence of solutions of the outer and inner problems 

In this section, we show that the solutions Ui and V^ of the previous equations exist 
under assumption (H3) on the spectral stability of the viscous shock profile. We 
first remark that the leading-order outer function u is exactly the solution of ^ 
which we want to approximate. Therefore, we first prove the existence of the V^ , 
and then we prove the existence of ui and all the Vl- Simultaneously, we prove the 
existence of the 5q. Similarly, we prove the existence of M2, V^ and 5{. 
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Construction of the approximate solution 



2.2.1 Construction at order 

In this section, we deal with the existence of u, which satisfy equations ©-(111) 
and matching condition ( |T0|) . The existence of u is exactly assumption (H2). Since 
M is a distributional solution, u verifies Rankine-Hugoniot conditions at each shock 
j- 

This property ensures the existence of the viscous shock profile which verifies 

Vl^-{f{V^)-X'^{t)V\ = Q (13) 
and the asymptotic conditions: 

For the existence for all t of such a profile and its properties, we refer to N. Kopell 
and L. N. Howard [U]. In this section, we just recall the convergence rate of the 
profile and its derivatives as — ±oo. Since u^~^ and are hyperbolic rest points 
for the ordinary differential equation (JT3l) . we have for some a; > and for any 

a e N, 

\d:V^{U)-dru^Ht)\ < e--l«l, V^GM, (14) 
\d^V^{^,t)\ < e^'^l^l, V^gM. (15) 



2.2.2 Construction at order 1: existence of and ui 

In this section, we prove the existence of ui, V7, (Oj ^*) such that 

ui,t + {df{u) ■ ui)x - dg{u) ■ ui = u^^, (16) 
V^^^ - {{df{V^) - X;) • Vl)^ = Vl + Vlbl, - g{V^), (17) 
Vi{tt) = u^{t) + uiHm - m) + 0(1), e ^ ±oo. (18) 

We first remark that these equations are linear. As in [9j, it is convenient to deal 
with bounded solutions. Therefore, we write 

u( = vi - d{ 

where D-[ is a smooth function such that: 



iuir{t) ife<-i, 
iui+{t) ife>l. 



Consequently, Ul solves: 

U^^ - {{df{V^) - Xj) ■ Ui\ = 5l,Vl + t), (19) 
f//(±oo,t)=^f(t)-<5^(t)nf(t) (20) 
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with 

From estimates ([ED, (|T^ . we deduce that h satisfies: 

hKi.t) = f^n^Ht) + {df{u^^)-X'^)uiHt)-9{u^^) + 0{e-^\^\),a> 0. 

And, since m is a smooth solution of (El), we have h e L^(R) and: h^{^, t) = 0(e-"l«l)- 
Integrating ( |T9l) with respect to ^ yields 

- (dfiVn - X'^)Ui = 5l,,V^ + t h\r^, t) dr^ + C^{t) (21) 

^0 

where C^{t) is a constant, only depending on t. 

Let us solve the problem (|T6l) -( l2T]) with matching condition (1201) . Following [9], 
we construct the solution of this system in two steps. First, for all j, we fix t and 
5o, we find U{ solution of (1211) with finite limits at ±oo. Since these limits are 
explicit and only depends on t, ^g, and C-' , we use the matching condition fpUj) to 
rewrite f lTUj) as a hyperbolic boundary value problem where Ui and 5q are the only 
unknowns. After solving this system, we use the previous construction to obtain Ul 
solution of ( 12T|) with matching conditions ( 120|) . 

So, we fix t and 5q for all j. With exactly the same arguments as in [9], we prove 
the existence of Ul for all j. Hence, using assumption (H3) on the viscous shock 
profile and theory of Fredholm operators, we show that U{ exists and the limits 
satisfy: 

hin Uiii^t) = -(d/(w^±) - X'^)-\5l,u^^ + + C^) 

where = J^°° h{ri, t) dt]. 

We now use matching conditions ( ISO]) to eliminate in these relations. Indeed, 
we have 

(d/(«^-+) - x;)(«i+ - 6^:^) = -{5i,u^^ + + c^-), 

and their difference is 

A^+u{^ - A^'u{- + - = ^ii^'^ui^ - A'-ui-) - {W+ - W-) (22) 

where A^^ = df{u^^) - X'-{t). 

Now, we have to solve ( !TB|) . (1^ . In order to find a solution of this system, we 
rewrite it by fixing the shocks. Since the shocks do not interact, we can define a 
change of variable Z (see Figure [2]) which is bijective, continuous in (x,t), piecewise 
linear in x and piecewise smooth: 

X — Xit^ L L 
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Xi X2 Xi + L 



Figure 2: Example of the change of variable Z in the case m = 4, for some t. 



We also define vi by 

Ui{x,t) = vi{Z{t,x),t). 
It follows from these definitions that vi solves 

vu + {Z^ df{u) + Zt){Z-\t, z),t)vu + Hz, t)vi - h{z, t) = (23) 

where z x = Z^^(t, z) is the inverse of x ^ z = Z{t, x), and, h and h only depend 
on Z, u, Uxi and u^x- 

We now use the fact that d/(-u) is diagonalizable, df{u) = P{u)^^D{u)P{u) so 

{Zx df{u) + Zt){Z~\t, z),t) = P{z, t)-'D{z, t)P{z, t). 

Since zeroth order terms do not play any role in the wellposedness issue, we consider 
the simplified system 

wu + Diz,t)wix = Hz,t), (24) 
^i+(pi+)-i^+ _ A^-(pj-)-^uj(- + 4(Mi+ _ u^-) = V{t) (25) 

with k and P known functions. Therefore, we have to solve this system on [0; L\ 
under periodic boundary conditions. Since and L correspond to the same shock, 
the periodic boundary conditions are in fact the shock conditions (1251) for j = 1. 

Equation (!^ is a hnear transport equation on Wu^i = 1, . . . ,n. Since, gener- 
ically, the existence of u smooth on [0; T*] ensures that the characteristics do not 
intersect on [0; T*], they can be used to build wi^i smooth between shocks, using the 
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initial condition. So, it suffices to verify that the conditions ( 125|) at the shocks are 
well-posed. We must therefore count the incoming and outgoing information at the 
shock. As we can see in Figure El for i < k, the incoming characteristics come from 
the right, so we obtain the value of . For i > k, the incoming characteristics 
come from the left, so we get w{~ . And for i = k, the sign of the eigenvalue change 
between two shocks: negative on the right of a shock and positive on the left, so, 
using again characteristic construction, wik is defined on the whole interval delim- 
ited by the shocks: we obtain ty-J^ and w^^.By this method we have built Wu on the 
right or left side of each shock. We now use the boundary conditions ( l25ll to obtain 
all the components of w{, and ^ot- 




Figure 3: Characteristic curves between two shocks, example with m = 2. We also 
plot Si and S2 which are sonic points for k-th eigenvalue, that means Xk{sj) = 0. 



Indeed, if we note by rj the i-th eigenvector of df{u), and Wi = PY^^aiTi, then 
a^"*" is known for i <k and for i > A; by our construction and we rewrite f l25p as 



or equivalently 



3+ j- 



V{t). 



1=1 
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This yields the hnear system on the unknowns a^"*" for i > for i < k, and 6^^ 

i>k i<k 

i<k i>k 

These equations have unique solutions if and only if the system obtained is invertible 
for all t, that is the Majda-Liu condition: 

Vj = 1, . . . , m, det(r^", . . . , r^I^, - m^", r^+^, . . . , r^+) 7^ 0. 

Using ^12j, our assumption (H3) implies Majda-Liu condition. To finish the con- 
struction of the approximate solution, we use again the characteristics. By this way, 
we have built a solution on the whole space R. To ensure the regularity of the so- 
lution far away from the shocks, we only need suitable compatibility conditions on 
the initial data. 

Finally, we have proved the existence of Ui and 6^^ for all j and for < t < T*. 
The previous construction give us for all j. We can then apply the same method 
to obtain the existence of V2 , Sl^ and U2, since the linear system has the same terms 
of maximal order. 

Remark 1. Since the construction of viscous shock profile only depends on the 
shock, we can use the previous construction even if u is not periodic. However, we 
will see in the following that the periodicity of u allows us first to obtain bounds on 
Ml and secondly to build the Green's function in Section [3l 



2.3 Construction of the approximate solution 

We complete the construction of an approximate solution of equation ([T]). First, we 
define a smooth function /i such that: 

= / if > 2, 
\ 1 if |x| < 1. 

Then, the approximate solution u^^p is defined as 



and ui„„ verifies 
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where t) = Yl^=i 9'f (^5 ^) is an error term given by 

qlix, t)={l- f,^) [{fm -f{u)-e df{u) -u.-e' df{u)-u, - ^ d'f{u) ■ K, m))^ 

- - g{u) - e dg{u) ■ Ui - dg{u) ■ ^2 - y d'^g{u) ■ 

ql{x,t) =^^^[(f(P^)-f(V^)-edf{V^) ■ Vi -e'df{V^)-Vi-'-^d'fiVn-iVi,Vn). 
-{g{P^)-g{V^)-edg{V^)-Vi) 

t) =f^l{P^ - on - sf^Ul^^ - on - 2e^,i{P^ - 0% + - /(O^)) 

+ /(^^P^ + (1 - f^nO% - if^'fiin + (1 - f^')fm). 

- gifi^^^ + (1 - fi^on - if^'gin + (1 - f^')9m), 

and /i-^ = /i j . 

We now want to prove that u^pp is a good approximation of u^, that means 
— M^pp —J- when e — )■ 0. Therefore, we define = — w^^p which solves: 

+ (d/«pp) ■ w)^ - ew^^ - dg{ulpp) ■ w = -q' + Qi«pp, w) - Q2«pp, w)^ (26) 

with 

QiiKpp, w) = g{w + <pp) - ^«pp) - d^«pp) ■ w 

and 

<52«pp, w) = f{w + <pp) - /«pp) - d/«pp) ■ w 
which are at least quadratic terms in w. 

2.4 Estimates on the error term 

To end with the proof of Theorem [21 it remains to compute the estimates on the 
error term q^ . As in |3j, we can estimate the support of functions g|: 

supp(gD C {(x,t) : \x-X,{t)\>e^], 
supp(g|) C {(x,t) : \x - X,{t)\ < 2e^}, 
supp(g^) C {(x,t) : < \x-Xj{t)\ < 2e^}. 

To obtain estimates on gf and their derivatives, we first need to fix the shocks by 
a smooth change of variable. This manipulation cannot be avoided for the estimates 
on qlt, qltt. So, we define 

m , L \ 

^{z,t) = z + 5^a,(^,t) X,(t) - (j - 1)- - eb^(^) (27) 
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where aj{-,t) are smooth functions, such that Xlj % = 1, is increasing, ip^ > 
and aj{-, t) = 1 on a neighbourhood [(j — 1)^ — r; (j — 1)^ + r] of (j — 1)^- We recall 
that in assumption (H2) we have supposed that \Xj^i — Xj\ > 2r, which ensures 
the existence of such a if. 
With the notations 

w{z, t) = w{ip{z, t),t),ulpp{z, t) = ulpp{ip{z, t),t), q'{z, t) = q'{(p{z, t),t), 

equation (l26l) becomes: 

wt + —{df{ul ) ■ w), - (— + ^ ) w,-e i^wA - dg{ul ) ■ w 

= -f + Qi{ulpp,w) Q2{ulpp,w),. (28) 

We now prove estimates ([n])-©"®-®- Using the fact that / and g are smooth 
and Ui is piecewise smooth, with discontinuities only at x = Xj{t), we have the 
following estimates for gf: 

llgilU- < 

where C is a constant which does not depend on e. Integrating this inequality, we 
get: 

||q^iIIli(0;L) < Ce'^. 
Moreover, if e is small enough (such that 2e"' < r), 

&(z,t) = (1 - . . 1. + e^-'Si^^ ( 



where the [. . . ] is dominated by e^. Therefore we have 

||^iJl°° < C^^, ||^iJli{0;L) < Ce^. 

Similarly, we prove the same estimates for ql^^. Nevertheless we cannot have such 
estimates for gf^. Indeed, we have 

fc(M)-(i-A')i...i.-.-v( ^''-";-^^'" 

so we just have 
Similarly, we prove 

||gLtlU°° < Ce^~^, \\qut\\L^O;L) < Ce^, 
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Estimates on the error term 



<Ce'-'\\\qU\LHo;L)<Ce 
Then, we compute estimates for 

=0 < Ce^ 
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LzIUmO;^) 

Eventually, we use matching conditions ( ITU|) - (|TTI) - (fT^ to prove the estimates on 
g|. Indeed, the properties of viscous shock profiles provide that terms o(l) can be 
replaced by e~"'^' in the matching conditions with a a positive number. So, we have 
for 2 — (j — 1)— < r: 



^-(i-l); 



t] +eVi 



e J - I ■ - 2 \ e 

L 



,t 



''-U-i)jr.-^si 



+0 e 



~^ 2^ XX i^) 



1^-0-1)^1 



and, using Taylor expansion for Ui{Xj + + (j — 1)^ — eSD), 

0^{if{z,t),t) = u{z + X,-{j-l)j^-e6^,t) 

+eui{z + X,-{j-l)j^-e6^,t) 
+e^U2{z + Xj-{j-l)j^-e6^,t) 

= +u^Ht) + uiHt){z - U - 1)^ + e6i + eH{] 



+eu{^{z-{j-l)j^-e6i)+e^u 
+0{e^ + {z-{3-l)i-e5lf) 



Since 



we have 



supp(g3) C <^ {z,t) : < 



z-{j-l)--e5^ 
m 



< 2e^ 



0\e'+(z-{j-l)--e5l 

m 



0{e 



37 ^ 



and we obtain the estimates: 

\mL-AmL-,mu\\L^<ce'\ 
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\\q3\\L^{0;L), Msth^iOiL), Mmh^iO-.L) < Ce^'^ , 

||gLJU°° ^ II5'LJIli{o;l) < C'^'^- 

This ends the proof of Theorem O 

3 Estimates on the Green's function 

We now consider the hnear operator 
Uw = wt + — (df{ul ) -ipt + ) ■ - e\w^^ 

+ ( ( ^ (^dfiulpp) -ipt+ ) - d^(Mapp) ) ■ 



The aim of this section is to prove the following theorem: 

Theorem 4. There exists a Green's function G'^{t^ r, z, y) of the linear operator U 
defined for 0<T,t<T*,z,y^W such that G'^{t, r, z,y) = if t > t and 

sup [ [ \G%t,T,z,y)\dzdt + y/e sup [ [ \d,G'{t,T, z,y)\dz dt < G 

y,T<T'- Jo Jr y,r<T* Jo Jm. 

(29) 

where G is positive and does not depend on e. 

To find estimates on this Green's function, we use approximations of the Green's 
function both near the shock and far away from the shocks. 

First, we recall the method of iterative construction of the Green's function of 
E. Grenier and F. Rousset Then, to approximate the Green's function near the 
shocks, we recall the result of K. Zumbrun and P. Howard [TT] about Green's function 
for pure viscous profile problem. Far away from the shocks, we use characteristic 
curves to build some approximate Green's functions. Finally, we combine all these 
Green's functions to obtain an approximate Green's function of L"^ and we find 
bounds on the error terms. 

3.1 Method 

Here, we recall the method used by E. Grenier and F. Rousset in J4j. We want to 
construct an approximate Green's function G^^^ of U in the form 

N 



G'appit, ^, z,y) = Y^ Sk{t, r, z, y)Hfc(r, y), 



k=l 
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where Sk are Green's kernels which satisfy (EO]) and Hfc e C°°([0, T*] x M, £(M")) are 
such that 

||nfc(t,x)t;|| < Va; > 0,t G [0;T*],t; G 

and 

N 

J]nfc = id. 

fe=i 

We next define the error Rk{-,T, -ji/) = L^Sk{j,y) for k = 1,...,N, and the 
matrix of errors: M{Ti,T2) = {(Tki{Ti,T2))i<k,i<N with 

akiiTi,T2) = sup [ [ \Uk{t,z)Ri{t,T,z,y)\dzdt. 

Ti<r<T2,J/Gsuppni JTi Jm. 

Thanks to Theorem 2.2 of [4j, we just have to prove that there exists 62 > such 
that < T2 — Ti < £2 iniphes 

hm MP{Ti,T2) = 0. 

p— >-oo 

Remark 2. Since we consider the error in a matrix, we need to consider a finite 
number of Green's kernel Sk- We will see in the following that the number of these 
kernels is proportional to the number of shocks per period. So, this method does 
not allow us to treat the case of a non-periodic perturbation. 

3.2 Near the shocks 

Near the shock j, we can approximate 

V ^z 



by 



df I ( - — t\ \ -X' 



3 



and we forget zeroth order term. Therefore, we search the Green's functions for the 
linear operators 

which depend on j and t < T*. As in [9], we remark that these Green's functions 
Gf^it,z,y) verify 
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where G^^^ is the Green's function related to the operator 

L> = dtw + {dfiV^iz, r)) - X;{t))w, - w,,. 

In [TT], K. Zumbrun and P. Howard obtained estimates on the Green's functions 
which will be useful to obtain estimates for our operators. Let us denote by a|^(r), 
and r^^(r) the eigenvalues and the associated eigenvectors of d/(M-'^(r)) — Xj(r). 

Proposition 1. Under hypothesis (H3), we have 



y m mj J-^ I 




r\ (r)x^<o 



(30) 



a.Gf ft,^ + (j-l)-,y + (j-l)-')= 

\ m mJ ^-^ \ 



j,a^+(r)>0 \ 




rl {t)xz<q 



(31) 

where M and o are positive constants, and x designs characteristic function. More- 
over, O 's are at least linear forms, locally hounded in y and uniformly hounded in 
it,z). 

3.3 Far away from the shocks 

As in the previous section, we do not search Green's function for but for the 
approximate operator defined by 

= wt + — [df{u{if{z,t),t)) - (fit + e'^] ■ - ^w^^. 
^z\ ^IJ 
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We remark that, as in the previous section, we forget the terms in w. 
Recall that df{u{x, t)) is diagonahzable for all x, t, so we can write 

d/K(^(z, t),t)) = P{u{^{z, t),t))D{u{if{z, t),t))P{u{if{z, t),t)) 



-1 



with D{u{(p{z,t),t)) = diag{Xi{u{(p{z,t),t))). 

To obtain approximation of the Green's function between two shocks, we define 
j approximate problems on IR, with continuous solutions. First, we set 



Xi{u{X+{t),t)) iize 
XMXi+iit),t)) iize 

Then, we want to find approximate Green's functions for the scalar operators 



To do so, we define characteristic curves r, y) by 

dtxl (t, ^, y) = K ixi {t, r,y),t), t>T, 
xl{r,r,y) = y, 

and the approximate Green's functions 

Gl{t,r,z,y) = ^M^e.p ( i^M ' ^^i^^-Mv.rW 



v/47r£(i - r) \^ ^e{t - r) 

We easily compute the error committed here 

Ij>G{ = {Xliviz, t),t)- Xiijiit, r, ^{y, T)),t))Gi(t, r, z, y). 

Before we build the whole Green's function, we introduce some notations. First, 
we write 

G^=diag(G^). 

In the sequel, we need to distinguish at each shock the outgoing waves to the in- 
coming waves. We define 

£)-jn _ (iiag(0, . . . , 0, 1, . . . , 1), with A; — 1 unit coefficients, 

jj-out _ cliag(ij . . . , 1, 0, . . . , 0), with k — 1 null coefficients, 

jj+tn _ fiiag(l^ . . . , 1, 0, . . . , 0), with k unit coefficients, 

jj+out ^ diag(0, . . . , 0, 1, . . . , 1), with k null coefficients, 

so that + D^°^^ = Id. 

Finally, we define the projections 

p±^«(t,^) = P{t,z)D^'''P{t,z)-\ 
V^'^\t,z) = P{t,z)D^'"'^P{t,z)-\ 
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3.4 Approximate Green's function 

Since the shocks are non-characteristic Lax shocks, we have the following inequality 
on a neighbourhood of each shock j = 1, . . . ,m: 

\\i(u(z,t)) - X.j(t)\ > C > if (j - 1)— -4:r]<z<(j - 1)— + 4r7, i = 1, . . . , n. 

m m 

We can assume that t] is such that ir] < r so that aj = 1 in flTTl) on the previous 
neighbourhood. 

Furthermore, we need some cut-off smooth functions 

+ / if 2;< 1, , X / 1 if ^< -2, 
^ ^^) = \ lif.>2 ^^^^ (") = \ Oif.>-l. 

We also assume the cut-off function fi already used to read as = {1 — K^){1 — K^). 
We can now build an approximate Green's function in the form 

m 7 
j=l k=0 

where the Green's kernels are periodic with period (0, 0, L, L): 

slit, r, z,y) = J2 slit, T,z + lL,y + IL) 

with 

~Sl{t,r,z,y) = /. (^ "~^\~'^- ^ (^^^37^) ^r'^t~r,z,y), 
Sl,it, T, z,y) = fi (^^^ l^i^^Zlli j p(t, ^)Z}+-*G^(t, r, z, y)P{T, y)-\ 
~Si{t,T,z,y) = /X (^^^ ir+ [ ^L'e^ l Pit,^)D^'''G^it,T,z,y)PiT,y) 



-1 



Si{t,T,z,y) = fi(^^^^^ K- (^^^Pj P{t,z)D-^-G^{t,T,z,y)P{T,y)-\ 
Syt, r, z,y) = fi (^^^ K- i j P(t, r, ^, |/)P(r, 
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and the projectors are also periodic: 

ni{T,y) = J2Kir,y + iL) 

with 



77 / \ M2£ 

y \^ M2£ \^ M3£ 

ni(T,i/) = K+ ( ^~^^~^^^ ^ + K- (^^^^j - 1> 

It appears that all the Green's kernels can be written in the following form: 

Si{t,T,z,y) = T{z)S{t,T,z,y) (32) 

where T is a truncation function. 

We will choose the three constants Mi,M2,M3 at the end of the estimates on 
the error matrix so that they verify 

4Mi < M2 < ^Mg. 

These inequalities are necessary to have 
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and 

Under these notations, describes the viscous dynamic at the shock j, S{ the 
creation of outgoing waves in a vicinity at the right of the shock j, 5*2 the creation 
and propagation of outgoing waves away from the shock j, at its right, 5*3 the 
creation and propagation of incoming waves at the right of the shock j. S\ describes 
the propagation of the waves between the shocks j and j ' + 1 Moreover, the kernels 
5*7, S'g, Si are the symmetric of respectively S{/S2, for the left of the shock j + 1. 
We summarize this splitting in Figure HI 




Figure 4: Summarize of the splitting by Green's kernels. 



3.5 Bounds on the error matrix 

As said in Subsection 13.11 to prove Theorem HI it remains to prove that Ai"^ con- 
verges to when p goes to oo. Since the coefficients of M. are non-negative, it 
suffices to prove that M. is bounded above by an other matrix which has the "good" 
convergence. 

To bound the error terms, we use the same method than F. Rousset in [0]. We 
split all the error terms into two parts: the truncation of the error on the kernel and 
the commutator: 

where, with the notation of fl32p . 

El{t, r, y) = T{z)L'S{t, r, z, y) and r, y) = [L^ T{z)]S{t, r, z, y). 

Lemma 1. We have the estimates 
at shock j : 

\\'^\y-U-l)L/m\<2M2eEQi{t,T, Z,y)\\Lo<>^^Li^^ < Ci(T + e) , 
\\'^\y-ij-l)L/m\<2M2e'^±z>oV^utEl2i'^^'^^ ^^y)\\L^y,Ll^ < ^2, 

\\hyHj-'^)LM<2Ahe'^±z>oV^El2it,T,z,y)\\Lr^^^^^^ < C3 + C2T ; 
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for the outgoing waves: R{, R2, R^, Ri'i ■ Let M > M2. we have: 

II hj-{J--^)L/m>MeEn {t, T, Z, y) , , II lj/>M£^21 ^> y) , 

<C4(r + £27-l) + C'5, 

II ly^(j-i)L/m<-MeEii (t, T,z,y) ^ , \\'^y<-MeEn 'r,z,y) Wl^^^l] ^ 
Uy-{j~i)L/m>MsEi2{t,r,z,y)\\L^^^Lj^^, Uy>MeEi2{t,r, z,y)\\L^^^Ll^<C5+Ci{T+e) , 

II h-U-i)L/m<~MeEl2 (t,T,Z,y) || L«>^,Li ^ , || ly<-MeE72 (t,T,Z,y) Wl^^^l] ^ 

<C5 + Ci(r + £); 

for the incoming waves: i?3, R^ 

II ^y-{j-l)L/m>AIeEii {t, T,Z,y) ^, II ^y<-MsEl^ (t, T,Z,y) ^ 

<'C6(T + £27-l)+C7, 
\\'^y-ij-l)L/m>MeEi2{t,T, Z,y)\\Lo.,^^Li ^, || ly<_Af£-E^2 ('^J ^) ^5 Z/) II L?°j,,LL 

' <Cs + C^iT + e); 

between two shocks: R\ 

\m\Lr.,,Ll.^<C,{T + e), 

where Ci is locally bounded in M2, M^,C2 is locally hounded in M2 uniformly in M3, 
C3 depends only on M2 and M3 and goes to as M3 — )■ +00, C4 is independent of 
Mi,M2 and M3, C5 depends only on M and goes to as M ^ +00, Cq is locally 
bounded in Mi, C-j goes to as Mi — > +00, and Cg is bounded uniformly in Mi. 

Proof. We do not give here the complete proof of the lemma. Mainly, it deals with 
terms that are not treated in [0]: zeroth order terms and terms related to (f or rj. 

First, we consider the error at the shock j. On the support of i?oi, we have that 
(Pz{-,t) = 1 and (pt{-,t) = Xj(t) — 661 (t) so we obtain 



2ri r \ MsE 



hG',\t-T,z,y) 



+ {dfiKJ - df ^V^ (^^Jl_i)M, ^ + x;iT) - X'^it) + e5^^ Gli 

where /i is a bounded function of t, z. Hence, using the fact that d/ and Xj are 
smooth and the expression of w^pp, we have 

li^^il < (^-^^^^^^ ^ (^^^37^) [(l^--l +^)l^'^l + l^'^O 

where C is locally bounded in M3. The calculations of [S] give directly a bound for 
the first term, in G^i- So, it only remains to bound the integral: 

/ / \Gl\t-T,z,y)\^z^t, 

Jt J {j-l)L/m-2Mze 
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and, thanks to Proposition [H to estimate: 

.M,. 1 / (^-a(t-r)r- \ 

J-Ahe Vs(t - t) V Me{t-T) J 

Using the classical change of variable z' = ^7"^*^""^ , we obtain the bound 

y/Me{t-T) 

j-T r{j~l)L/m+2Mse 

/ / \G'^^{t-T,z,y)\dzdt<CT 

Jt J {j-l)L/m-2Mse 

SO this concludes the estimate for Eq-^^. For Eq2, the linear term in Gf^ disappears 

in the commutator. Hence, there is no change with the proof of F. Rousset. 

For the estimates on Rl,i = 1, 2, 3, 5, 6, and 7, we first remark that v?^ = 1 on the 

support of the errors, so E^^ is bounded as in jH], except for the zeroth order term 

which is treated as in the case of i^gi- So, we only consider the estimate on i^/j- We 

first remark that the support of this error is not of size e. Indeed, the truncation 
Z^— (7— 1) — \ 

/X ( ^ — —j adds some error terms with support of size Arj. Though, we have to 

bound: 



-e -^Li"K+PD+°'''G^p-^ + —Li'K^{PD+°'''&),p-^ 
Ari^ 2r] 



dzdt. 



The two terms in G^ are bounded by GT, and the term in {PD^°'^^G^)z is bounded 
by G{e + T). Indeed, for 1 < i < n. 



T /•(i-l)|r+4r, 



\Gl\ 



diiciuic ,0 



Hence, using the change of variable z' = xAj^^V''^)) ^ obtain 

y/4e{t-T) 



I I 



G'\<G I I e""dz'dt < GT. 



Similarly, we have 



r J[j-i)i.+M,e ^ATxe{t - t) \ 4e(t -r) \ 2e(t - r 




exp -^.^^ 1 dzdt. 
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And with the same change of variable, we obtain: 

pT .0-1)^+4,, r'^ 1 r ,2 

e / / \GU<C^ / ^= / \z'\e-' Az' At < C{T + e), 

which gives the estimate for £'^2- 

Since there is not new difficulty in the proof of the estimates on we do not 
develop it here. □ 



We now use Lemma [T] to bound the matrix M.. Since the Green's kernel depends 
on the shock, we note 



T2 



al{T^,T2)= sup / j \Yi\{t,z)Rl{t,T,z,y)\dzdt. 

Ti<r<r2,j/6suppn^ -JTi 



Since two shocks do not interact, the error coefficients a^^i vanish for i ^ j and kl ^ 
and for \i — j\ > 1. So the error matrix Ai is bounded: 



/ Ml M2 
M3 Ml M2 



M < 



\ M2 



M3 \ 



■•• M2 

M3 Ml / 



where M2 is null except on the first column, and M3 is null except on the first line. 
Moreover, using Lemma [Tj and it was done in [9], we can choose a < 1/2 and Mi, M2 
such that when M3 — )■ +00, e, T — )■ 0, the matrices tend to 



Ml ^ 



a 


C 






\ 


a 


a 








C a 


a 








a 


a 








a 


C 


C 


a 








a 


a 








a 


a 








a 


a 


) 



( 



M2 ^ 



c 
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C 



a \ 



V 



So that Ai ^ M. where the eigenvalues of M. are and 2a. We can now conclude 
that we have Ai^ — ?■ when p — )■ oo. This ends the proof of Theorem |H 



4 Convergence 

The purpose of this section is to prove Theorem [3] and to conclude the proof of 
Theorem [H For this sake, it remains to show that solution of 



w{z,0) = 



w 



1 



(33) 
(34) 



vanishes as e — j- 0. In the previous section, we obtain estimate (!29l) on the Green's 
function of operator U . We recall that Qi and Q2 are at least quadratic in w and 
(f verifies inequalities ([6]), (JTj), (jH]), and ([9]). 

As in |1] and we use standard arguments for parabolic equations. First, 
we remark that local existence of a smooth solution w for (!HH|) - flM|) is classical. Then 
we define 



= sup{Ti G [0; T*], 3w solution on M x [0; Ti), EiTi) < 1}, 



where 



^(Ti 



^0 



\w\ 



+ 



+ 



37-a-l/2 ' g-37-2a--l/2 



£37- 



2a-l 



+ 



+ 



g-37-2a-l g-37-3a-3/2 ^37 



-2a-2 



with a > and 7 G (2/3, 1), chosen later. 

Before estimating the L^— norms, we define the notation: 



\\w 



|l = ||w7||li{(0;T'=)x[0,L)), II^IIoo = II^H loo((0;Ts:)x [0,L)) 



dzdt 
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First, we use Fourier coefficients c„ to estimate || Indeed, 
w{z,t) = [ wt{z,s)ds= [ ^c„(u't(s))e*"^ds 

= / co{wt{s))ds+ / e ds 

"^0 -^0 nez\{o} ^ 

\\w\\oo < C{\\Wth + Wi^zztWl) < C£-^^-3a-3/2 

when e < 1. Hence, llt^lloo tends to as e goes to if 7 and a are such that 
7 - a > 1/2. 

We deduce that to prove that the time existence is T* and that we have the 
convergence, it remains to prove that = T* for e small enough. In the sequel, we 
suppose that < T* so E{T^) = 1. 

Then, we can use the estimate on ||w||oo to bound w and its derivatives. First, 
using we have 

w{t,z)= I I G%t,T,z,y) (-q' + Qi{ulpp,w) - —Q2{ulpp,w)zj {T,y)dydT. 

Also, we have for periodic function tl) the estimate 

^ G\t, r, z, y)d^ ^{y) dy = O (^j^ G%t, r, z, y)dz^ ^(y) dy^ . 

Therefore, we deduce 

llu-lli < Ce^^ + CII^DlUll^lli + CII^Dlloodl^lli + ll^.lli) 

so 



F 1 



This can be made smaller than 1 as e — )• if a and 7 are such that 87 — 3a — 2 > 0. 
We now take the z derivative w and obtain an expression of Wz'- 

Wz{t,z)= I I Gl{t,T,z,y) (-f + Qi{ulpp,w) - —Q2{ulpp,w)A {T,y)dydT 



and the estimate: 

ll^^zlll 



e 



37— a— 1/ 



Differentiating equation ( l33l) with respect to t or 2;, we obtain equations verified 

by Wt and Wz- 

L'wt = {-q" + Qi{ulpp,w) - Q2iulpp,w)z)^ + h{w,wz,ewzz), 
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where /i and I2 are continuous linear forms, uniformly bounded with respect to e. 
Thus, using again the Green's function, we get the inequalities: 



^87-20-1/2 



< C(£2-+V2 + ^37-3^-3/2 (^a+1/2 + 1 + ^" + ^"1/2) + ^«+l/2 ^ ^ ^1/2 



< C(£2"+V2 + ^37-3a-3/2^^a+l/2 ^ ^-1/2) _^ ^ ^1/2)^ 
i^|k < C'(^2.+l/2-7 + ^37-3.-3/2(^^+1/2 ^ ^-1/2) ^ ^ ^l/2y 

In the last inequality, the right-hand side can be made smaller than 1 as e — )■ if a 
and 7 are such that 2q; + 1/2 — 7 > 0. 

Differentiating again equation ( |33l) with respect to t, z or we obtain the 
equalities: 

where 1^ and ^4 are continuous linear forms, uniformly bounded with respect to e. 
As seen before, we deduce the inequalities: 

ji-3l-3/2 ^ ^^(e^^+l-T + ^37-3a-3/2(^2a+l ^ ^-1/2) ^ ^2a+l ^ ^ ^1/2)^ 



^^^^lll < (^(^^2a+3/2-27 _^ ^37-3a-3/2 (^^a+3/2 _^ £--1/2) _)_ ^Q+3/2 _^ ^1/2^| 



-2a-2 



The both bounds can be made smaller than 1 as £ —t- if a and 7 are such that 
3a + 1 - 7 > and 2a + 3/2 - 27 > 0. 

We can now verify that there exist 7 and a checking all the previous conditions: 
they define a non-empty trapeze in the plane a, 7. Hence, we have proved that for 
e small enough, we have E{T'^) < with /3 > 0. Consequently, we can not have 
rpe ^ such an e. Moreover, using the change of variable (p, we return to 

w = — u^app- Thus, we have 

ihiloo < 

and 



Wu^-u" II < r<p37-3Q-3/2 
II " "applloo — 

Inequality E(T*) < 1 also gives 

Wwh^iL^) < \\wt\\i < C\\wth < 
28 



so 

lU,"^ II <r r'c-37-2a-l/2 

11^ - ^applU°°(Li) < <--^ 

This concludes the proof of Theorem O 

It only remains to prove Theorem [H Since H^app — — > 0, we have the 
convergence in L°°((0; T*), -L^(0; L)). Moreover, the fast convergences of the viscous 
shock profiles give the last point of the theorem. 

5 Conclusion and perspectives 

In this article, we have proved the persistence of solutions of the inviscid equation ([2]) 
close to roll-waves by adding full viscosity. One of the main assumptions that we 
have taken is the periodicity of the solution of (E]). This one is not necessary in 
the construction of the approximate solution w^pp, but it gives that U\ and U2 stay 
bounded (because periodic). An idea to weaken this assumption would be that 
solution uoi ^ approximates the roll- wave as |x| goes to infinity (in particular, the 
shock curves would be closer as goes to infinity). Moreover, the periodicity of u 
allows us to use the method of [1] to construct the Green's function of U . Indeed, 
in this step, the number of Green's functions that we consider is proportional to 
the number of shocks. In the periodic case, taking into account the repetitions, it 
returns to a finite number of periodic Green's kernels. Thus, we can write a matrix 
of errors, and deduce Theorem H] on the existence of the Green's function relative to 
U and estimates on this Green's function. Another way to hope to obtain a finite 
number of Green's functions could be to assume that u coincides with the roll-wave 
outside a bounded domain. 

Furthermore, Theorem [T] proved here is valid in the case of an artificial viscosity. 
Therefore, one should also study the persistence in the case of real viscosity as 
presented to the system of Saint Venant ([3]). 

Finally, one can also be interested in what happens in the multidimensional case. 
For this, we could build on work already done in the case of a single multidimensional 
shock, based on the study of Evans' functions at each shock [5]. 

Acknowledgment. The author wishes to thank Pascal Noble for suggesting the 
problem and for fruitful discussions. 
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